Refining new-physics searches in B — >• Dtv decay with lattice QCD 
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The semileptonic decay channel B — > Dtv is sensitive to the presence of a scalar current, such as 
that mediated by a charged-Higgs boson. Recently the BaBar experiment reported the first obser- 
vation of the exclusive semileptonic decay B Dt~v, finding an approximately 2a disagreement 
with the Standard-Model prediction for the ratio R(D) = BR(B -» Dtv)/BR(B -> Dlv), where 
£ — e,[i. We compute this ratio of branching fractions using hadronic form factors computed in 
unquenched lattice QCD and obtain R(D) = 0.316(12)(7), where the errors are statistical and to- 
tal systematic, respectively. This result is the first Standard-Model calculation of R{D) from ab 
initio full QCD. Its error is smaller than that of previous estimates, primarily due to the reduced 
uncertainty in the scalar form factor fo(o 2 ). Our determination of R{D) is approximately la higher 
than previous estimates and, thus, reduces the tension with experiment. We also compute R(D) 
in models with electrically charged scalar exchange, such as the type II two-Higgs doublet model. 
Once again, our result is consistent with, but approximately la higher than, previous estimates for 
phenomenologically relevant values of the scalar coupling in the type II model. As a byproduct 
of our calculation, we also present the Standard- Mo del prediction for the longitudinal polarization 
ratio P L (D) = 0.325(4)(3). 



PACS numbers: 13.20.He, 12.38.Gc 

Motivation. - The third generation of quarks and lep- 
tons may be particularly sensitive to new physics asso- 
ciated with electroweak symmetry breaking due to their 
larger masses. For example, in the minimal supersym- 
metric extension of the Standard Model, charged-Higgs 
contributions to tauonic B decays can be enhanced if 
tan (3 is large. Thus the semileptonic decay B — > Dtv is 
a promising new-physics search channel [1-13] . 

The BaBar experiment recently measured the ratios 
R(D^) = BR(B -» D^tu)/BR(B ->■ D^*Hv), where 
£ = e,[i, and reported excesses in both channels which, 
when combined, disagree with the Standard Model by 
3.4cr [14]. BaBar also interpreted these measurements in 



terms of the type- II two-Higgs-doublet model (2HDM II) 
and claimed to exclude the theory at 99.8% confidence- 
level. In this work, we update the prediction for R(D) in 
the Standard-Model and in new-physics theories with a 
scalar current (such as the 2HDM II) using unquenched 
lattice-QCD calculations of the B — > Div form factors 
foil 2 ) and f+(q 2 ) by the Fermilab Lattice and MILC col- 
laborations [15]. This is the first determination of R(D) 
from ab initio full QCD. 

With lepton helicity defined in the rest frame of the 
virtual W boson, the general expressions for the differ- 
ential rates for semileptonic B — > D£v decay are given 
by 
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where the subscript denotes helicity, and q = (pe + p v ) 
is the momentum carried by the charged lepton-neutrino 
pair. The total semileptonic width is the sum of the 
partial widths, r to t = (r+ + At tree-level of the 

Standard Model electroweak interaction, the scalar- and 
tensor-exchange couplings are Gs — Gt — 0, while the 
vector coupling is Gy — GpVij. In the infinite heavy- 
quark-mass limit, the form factors f + (q 2 ) and f2(q 2 ) are 
related via 

h{q 2 ) = -U(q 2 ) Ml ~ 2 Ml [U( q 2 ) Ml 2 )] ■ (3) 

The form factor f2(q 2 ) is only relevant for theories with 
tensor currents, however, which we do not consider here. 

Because the Standard-Model positive-helicity contri- 
bution to semileptonic B — > D decay is proportional 
to the lepton mass-squared, it can be neglected for the 
light lcptons £ — e, fi; thus experimental measurements of 
B — > D£v decays are sensitive only to the vector form fac- 
tor f + (q 2 ). On the other hand, the differential rate for 
B — > Dtv is sensitive also to the positive-helicity con- 
tribution and, hence, fo{q 2 )- Existing Standard-Model 
estimates of dT-/dq 2 for B — > Dtv have relied on the 
kinematic constraint /o(0) = /+(0), dispersive bounds 
on the shape [16], relations from heavy-quark symme- 
try, and quenched lattice QCD (neglecting u, d, and s 
quark loops) [17, 18]. See Refs. [10-13] for details. In 
this letter, we replace quenched QCD and heavy-quark 
estimates with a full, 2+1-flavor QCD calculation. In 
particular, we determine the following ratios within the 
Standard Model (where I = e,/x): 

R(D) = BR(B -> Dtv)/BR{B -> Dlv) , (4) 

P L {D) = (tB-^Dtv _ pB-vDr^J /Tf^ DTV . (5) 

These quantities enable particularly clean tests of the 
Standard Model and probes of new physics because the 
CKM matrix elements and many of the hadronic uncer- 
tainties cancel between the numerator and denominator. 

Lattice- QCD calculation. - Here we briefly summarize 
the lattice-QCD calculation of the B — > Dlv semileptonic 
form factors f+{q 2 ) and fo(q 2 ) [15]. Our calculation is 
based on a subset of the (2+1) -flavor ensembles generated 
by the MILC Collaboration [19]. We use two lattice spac- 
ings a w 0.12 and 0.09 fm, and two light-quark masses 



at each lattice spacing with (Goldstone) pion masses in 
the range 315-520 MeV; the specific numerical simulation 
parameters are given in Table I of Ref. [15]. 

This relatively small data set is sufficient for ratios 
such as those studied here and in Ref. [15], given the mild 
chiral and continuum extrapolations. We use the Fermi- 
lab action [20] for the heavy quarks (bottom and charm) 
and use the asqtad-improved staggered action [19] for the 
light valence and sea quarks (u,d,s). We minimize the 
systematic error due to contamination from radial exci- 
tations in 2-point and 3-point correlation functions by 
employing fits including their contributions, as described 
in Sec. Ill of Ref. [15]. We renormalize the lattice vec- 
tor current cj^b (and other heavy-heavy currents) using 
a mostly nonperturbative method [21] in which we de- 
termine the flavor-conserving normalizations nonpertur- 
batively. The remaining correction is close to unity and 
can be calculated in one-loop tadpole-improved lattice 
perturbation theory [22]. When extrapolating the lattice 
simulation results to the physical light-quark masses and 
the continuum limit, we carefully account for the lead- 
ing nonanalytic dependence on the light-quark masses 
at nonzero but small momentum transfer [23] includ- 
ing the effects of lattice artifacts (generic discretization 
errors and taste-symmetry breaking introduced by the 
staggered action) [24, 25]. The chiral-continuum extrap- 
olation results are plotted in the left panels of Figs. 6 
and 7 of Ref. [15]; for details see Eqs. (4.1) and (4.2) of 
the same work and the surrounding text. 

Figure 1 shows the results for f+(q 2 ) and fo{q 2 ) [15]. 
The simulated data are in the range w < 1.17 (to the 
left of the dashed vertical line), where w — (M§ + 
M|, — q 2 )j \2MbMd). In this region, we parameterize 
the w dependence of the form factors by a quadratic ex- 
pansion about w = 1, which works well for small re- 
coil. To extend the form-factor results beyond the sim- 
ulated recoil values (to the right of the dashed vertical 
line) we reparamcterize the form factors in terms of the 
variable z [27], and then extrapolate to large recoil us- 
ing a model-independent fit function based on general 
quantum-theory-principles of analyticity and crossing- 
symmetry. The functional forms used to extrapolate 
f+{q 2 ) and fo(q 2 ) are defined in Eqs. (5.1)-(5.6) of 
Ref. [15]. The fits are plotted in the left panel of Fig. 9, 
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FIG. 1. The form factors /+ (top, red) and / u (bottom, blue) 
from lattice QCD. The range of simulated recoil values is to 
the left of the vertical line. The filled colored bands show the 
interpolation/extrapolation of the numerical lattice data over 
the full kinematic range using the z parameterization. For 
comparison, the experimental measurement from BaBar [26] 
is shown as solid filled circles (using |V c i,| = 41.4 x 10 -3 [17, 
26]). 



and the results are given in the upper panel of Table V 
of the same work. As seen in Fig. 1, our result for f+(q 2 ) 
agrees very well with experimental measurements [26] 
over the full kinematic range. This nontrivial check gives 
confidence in the extrapolation of fo(q 2 ), which cannot 
be obtained experimentally and for which lattice-QCD 
input is crucial. In particular, lattice-QCD uncertainties 
are smallest near q 2 = (M B — M D ) 2 7 so the discussion be- 
low hinges principally on our calculation of fo(q 2 ) near 
this point, the validated /o(0) = /+(0), and a smooth 
connection between the two limits. 

We calculate the Standard- Model B — > Dlv partial de- 
cay rates into the three generations of leptons using these 
form factors and Eqs. (1) and (2) with Gs = Gt = 0, 
Gy = GipV* b . The resulting distributions are plotted in 
Fig. 2. To illustrate the role of the scalar form factor 
/o(<? 2 ), we also show the rates with only the contribu- 
tions from f+(q 2 )- Due to the significant helicity sup- 
pression, the differential decay rates into light leptons 
are well-approximated by a single contribution from the 
form factor f + (q 2 ). For B — » Dti/, however, the contri- 
bution from the scalar form factor fo(q 2 ) comprises half 
of the Standard-Model rate. 

Given the lattice-QCD determinations of f+(q 2 ) and 
fo(q 2 ) we can obtain the Standard- Model values for R{D) 
and Pl(D). These are the primary results of this let- 
ter, and we now discuss the sources of systematic uncer- 
tainty. In Ref. [15], many statistical and several system- 
atic errors cancelled approximately or exactly in the ratio 
/ B ^Dsi" /fB^Div studied tnere . Some of these do not 

cancel (as well) in R(D) and Pl(D), however, because 
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FIG. 2. Differential decay rates for B — >■ Dev (green), B — > 
D\w (blue), and B — \ Dtv (red) in the Standard Model. 
The black dot-dashed curves show the rates calculated with 

/o(g 2 ) = o. 

they affect f+(q 2 ) and fo(q 2 ) differently. 

Table I shows the error budgets for R{D) and Pl(D). 
The statistical error in R(D) is significant (3.7%) due 
to the different phase-space integrations in the numera- 
tor and denominator, whereas for Pl(D) the correlated 
statistical fluctuations largely cancel. For the same rea- 
son, the errors in R{D) arising from the extrapolation to 
the physical light-quark masses and the continuum limit 
(1.4%) and to the full q 2 range (1.5%), are much larger 
than for Pl(D). Wc estimate the error from the chiral- 
continuum extrapolation by comparing the results for fits 
with and without next-to-next-to-leading order analytic 
terms in the chiral expansion. We estimate the error from 
the z extrapolation by varying the range of synthetic data 
used in the z fit, including an additional pole in the fit 
function, and including higher powers of z. The specific 
chiral and z-fit variations considered are enumerated in 
Table VI of Rcf. [15] and discussed in detail in the sur- 
rounding text. The remaining sources of uncertainty in 
Table I do not contribute significantly to the quantities 
studied in Ref. [15], so we describe them in greater detail 
below. 

We determine the bare heavy-quark masses in our sim- 
ulations by tuning the parameters Kf, and k c in the heavy- 
quark action such that the kinetic masses of the pseu- 
doscalar B s and D s mesons match the experimentally- 
measured values [20] . In practice, it is easier to work with 
the form factors h±(w) on the lattice, which are linear 
combinations of /+,o(? 2 ) [15]. We study how the form 
factors h±{w) depend on Kf, iC by re-computing the form 
factors on some ensembles at values of Kb,c slightly above 
and below the default ones, and extracting the slopes 
with respect to Kf, jC . We use these slopes to correct our 
results for R(D) and Pl(D) slightly from the simulated 
k values to the physical ones, and conservatively take the 
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full size of the shift as the error due to K-tuning. 

We remove the leading taste-breaking light-quark dis- 
cretization errors in the form factors with the chiral- 
continuum extrapolation, and estimate the remaining 
discretization errors from the heavy-quark action with 
power counting [28] . We compute both the coefficient of 
the dimension 5 operator in the Fermilab action csw and 
the rotation parameter d\ for the heavy-quark fields at 
tree level in tadpole-improved lattice perturbation the- 
ory [20]. Then the leading heavy-quark errors in h + (l) 
are of O (a s (A/2m Q ) 2 ) and O ((A/2toq) 3 ), where A is 
a typical hadronic scale. Using the values a s — 0.3, 
A = 500 MeV, and m c — 1.2 GeV, we estimate that 
heavy-quark discretization errors in h + (l) are ^1-2%. 
At nonzero recoil, w > 1, there are corrections to h + (w) 
of O (a s A/2mg), but these are suppressed by (1 — w) 
because they vanish in the limit w — 1 by Luke's theo- 
rem [29]. We expect them to be largest at our highest 
recoil point w = 1.2, and estimate their size to be ~ 1%. 
Thus we estimate the uncertainty in h + (w) from heavy- 
quark discretization errors to be 2%, which leads to neg- 
ligible errors in R(D) and Pl(D). The leading heavy- 
quark error in h-(w) is of O (a s A/2mg), which we esti- 
mate with the input parameters above to be ~ 6%. To be 
conservative, we take the error in the ratio h_(w)/h+(w) 
to be 10%, which leads to small errors in R(D) and 
Pl(D). 

Our methods for computing B — > D transitions in- 
corporate the bulk of the matching of the lattice vector 
current to continuum automatically, leaving a factor p v » 
close to unity [30]. For R(D) and Pl(D), only the relative 
matching of the spatial and temporal components of the 
current matters, Pv\l Pv° ■ Although we have a one-loop 
calculation of p v o in hand, no nontrivial result for p V i is 

cb cb 

available. We take Pv\l ' Pv° = 1-0 i 0-2 to estimate the 
uncertainty from this source. In similar calculations, we 
have never seen p factors that differ from unity by more 
than 5%, so this range is extremely conservative. The 
uncertainty in the current renormalization factors leads 
to a small error in R(D), but is the second-largest source 
of error in Pl(D), after statistics. 

We also consider the systematic uncertainties from 
tuning the light-quark masses and determining the ab- 
solute lattice scale r\ , but these produce negligible errors 
in both R(D) and P L (D). 

Results and Conclusions. - We obtain the following de- 
terminations for the branching-fraction and longitudinal- 
polarization ratios for B —¥ D£v semileptonic decay: 



R{D) = 0.316(12)(7) . 
Pi(£>) = 0.325(4)(3), 



(6) 
(7) 



where the errors are statistical and total systematic, 
respectively. The value of R(D) is approximately la 
larger than the recent Standard-Model values obtained 
using estimates of fo(q 2 ) from Refs. [11, 13], but it is 



TABLE I. Error budgets for the branching fraction and lon- 
gitudinal polarization ratios discussed in the text. Errors are 
given as percentages. 



Source 


R(D) P L {D) 


Monte-Carlo statistics 


3.7 


1.2 


Chiral-continuum extrapolation 


1.4 


0.1 


^-expansion 


1.5 


0.1 


Heavy-quark mass (k) tuning 


0.7 


0.1 


Heavy-quark discretization 


0.2 


0.3 


Current Pyi /p v o 


0.4 


0.7 


total 


4.3% 


1.5% 



still 1.7cr lower than the recent BaBar measurement, 
R{D) = 0.440 ± 0.058 ± 0.042 [14]. The results for R(D) 
from Belle [31] agree with those of BaBar, but have larger 
uncertainties. Current experimental measurements of 
R(D) are statistics-limited, so the luminosities available 
at Belle II and SuperB should enable significant improve- 
ment on R(D) and possibly a determination of Pl(D). 

We also re-examine the interpretation of the BaBar 
measurement of R(D) as a constraint on the 2HDM II; 
the result is plotted in Fig. 3. For this theory, the scalar- 
exchange coupling in Eq. (2) is given by 



GpVcb- 



m e ( m c + m b tan 2 



(8) 



As soon as tan/3 > 4, m c in Eq. (8) can be neglected. 
On the other hand, evaluating Eq. (2) depends sensi- 
tively on m c /mb via the prefactor multiplying fo(q 2 ). 
We take m c /m b = 0.22, following Refs. [13, 14]. Our 
improved calculation of the scalar form factor fo(q 2 ) 
increases the prediction for R(D) by about la for all 
tan/3/M H + < 0.6 GeV" 1 . This sensitivity of R(D) to 
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FIG. 3. Lattice-QCD calculation (red) and experimental mea- 
surement [14] (blue) of R(D) vs. tm/3/M H ± in the 2HDM II. 
For the BaBar result, the dark- (light-) blue regions denote 
the lcr (2a) error bands. For comparison, the 2HDM II curve 
based on Refs. [11, 13] and used in Ref. [14] is shown in green. 
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differences in fo(q 2 ) suggests that one should be cautious 
in using indirect estimates of the form factors to con- 
strain new-physics models in other decay channels such 
as B — > D*tv. 

Inspection of the general formulas for the differential 
decay rates, Eqs. (1) and (2), shows that many new- 
physics explanations for the ~ 3ct tension between mea- 
surements of {R(D), R(D*)} and the Standard Model 
are possible. Even in the context of 2HDM, variants 
other than type II can have different phenomenology [32] . 
Lattice-QCD calculations of f+(q 2 ) and fo(q 2 ) can be 
used to provide reliable predictions for R(D) in any 
model given values for the couplings {Gs, Gy, Gt}- We 
note, however, that the Dalitz distribution for the lepton 
and D-meson energies in the i?-meson rest frame may be 
more sensitive to tensor interactions than R{D) [33]. 

The largest source of uncertainty in our determina- 
tions of R{D) and Pl(D) is statistical errors, which we 
expect to reduce with an analysis of our full data set in 
a future work [34]. The ratio R(D) is correlated with 
Pl(D) as well as other observables such as R(D*) or 
R(D S ) in many new-physics models; thus the pattern of 
experimental results for these quantities can help distin- 
guish between new-physics scenarios, such as those with 
and without a charged-Higgs boson [12]. We will present 
lattice-QCD results for R(D*) and P L {D*) in a future 
paper on the B — > D*lv form factors, and also note that 
we could easily obtain Standard-Model predictions for 
R(D S ) and Pl(D s ) if measurements of these quantities 
were possible with an T(5S*) run at a B factory. 

Given the present tensions between experimental mea- 
surements and Standard-Model predictions for both 
{R(D), R(D*)} and the leptonic branching fraction 
BR(B -> tu) [35-38], lattice-QCD calculations of B -> 
Dtv form factors and other hadronic weak matrix ele- 
ments can play a key role in revealing whatever theory 
beyond the Standard Model is realized in Nature. 
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